PRO DIANSP
C=[0,1,-1,-1,1,-1,1,1,-1,1,1,1,-1,1,-1,$
   1,1,1,-1,-1,-1,-1,-1,-1,1,1,-1,-1,1,$
   1,-1,-1,-1,-1,-1,-1,1,1,1,-1,1,-1,1,$
   1,1,-1,1,1,-1,1,-1,-1,1]

A=INTARR(53,53)
A=(C#C+1)/2
A[*,0]=0
A[0,1:52]=1

M=INTARR(95,95)
M[0:47,0:47]=ROTATE(A[0:47,0:47],2)
M[48:94,0:47]=ROTATE(A[1:47,0:47],7)
M[0:47,48:94]=ROTATE(A[0:47,1:47],5)
M[48:94,48:94]=A[1:47,1:47]

S=INTARR(189,189)
S[25,27]=18
S[65,70]=13
S[50,63]=16


B=RANDOMU(1,189,189,POISSON=1)


G=M+M-1
G[47,47]=1

D=CONVOL(S,M,CENTER=0,/EDGE_ZERO)+B

D0=INTARR(283,283)
D0[0:188,0:188]=D


S1=CONVOL(D0,G,/EDGE_ZERO,CENTER=0)

WINDOW,0,Xsize=283,Ysize=283
TVSCL,S1
ZOOM,FACT=2,XSIZE=566,YSIZE=566

END
